Reply to Comment on " Theory of Spinodal 
on"; Decomposition" 

S. B. Goryachev 

^ . Laboratoire de Metallurgie Physique, Bat. C6, 

r£^ I Universite de Lille I, 

59655 Villeneuve d'Ascq, France 
(N ■ 



(N 



Submitted in Physical Review Letters 06 July 1995 
>; : Revised 18 April 1996 



~^ ' In his Comment [|lj to my paper U Rutenberg notes that " the basis of my 

f^ . analysis of conserved scalar phase-ordering dynamics, to apply only the global 

^\ ' conservation constraint 

ON . / dx^ = N, (1) 

a is incorrect". Because "for physical conserved systems, which evolve by mass 

I transport, the stronger local conservation law embodied by continuity equation" 

^ ■ d^/dt = -Vj. (2) 

O' 

O ' This question needs clarification. 

Replying shortly, I say that my theory has no any contradiction with the 
- local conservation law (2). It describes the evolution of a^/ i/iree types of thermo- 

rS I dynamic systems: without conservation law (NCOP), with global conservation 

S ■ law (1) (GCOP) and with local one (2) (LCOP). In these all cases the order 

parameter (OP) Sip{x,t) = '0(x, t) — V'o(x) dynamics near the extremal V'o(x) 
is described by the same evolution equation 

d^j/dt ^ -rST/Sil; (3) 

with different thermodynamic potential functionals ^{ip}- For NCOP-system 
J-'{4'} = F{'4'} = /dx/l^/)}, where F{V'}is a free energy functional and /{i/j} 
is a specific free energy one. For GCOP- and LCOP-systems ^{1^} = ^{ip} = 
F — HgN, where the constant fj,g is a global chemical potential defined by (1) 
and equilibrium condition SJ^ / Sipl^^^^f^x) = 0. The difference between F and 
ri is that the conservation law (1) is taken into account. It should be noted 
that for GCOP and LCOP-systems (1) must be also used directly for selection 



of valid solutions of (3). So there is no difference between GCOP and LOOP 
evolution. That is why I have used the same notation (COP) for them in [pi. 

Writing evolution equation (3) I proceed from classical point of view that 
a system has a thermodynamic potential ^{V'} */ o^c' only if it is related to 
dynamical class of, so called, potential systems. For potential systems the OP 
evolution near the extremal V'o(x) is described by the equation (3). More gen- 
erally, ^{tp} is the system's Lyapunov functional for attractor -00, if .7^ < for 
tjj y^ ipo and T = ior tp = ipQ. Dynamics of the potential system consists 
of relaxation toward the minimum J^q = J^{iPq} and the Lyapunov functional 
existence guarantees a system's global asymptotic stability at ip = ipQ. In more 
physical language the condition T < ioi ip ^ ipQ and J- = for ijj = ipo is 
nothing more nor less than the second law of thermodynamics. 

Rutenberg states that (2) applies a supplementary constraint on the OP 
evolution. It is not right, because for LOOP we can always find from (2) and 
(3) the OP flux j , being a priori an unknown quantity, by writing a transport 
equation ^, |j] 

vj = rsn/Sip. (4) 

If rotj = 0, we can rewrite (4) as Q 

j == -rVM/(47r). (5) 

Here A4{x,t) = J dx'{[5n{tp{x',t)}/Sip]/\x — x'|} is a non-local characteristic 

of the system, so called, scalar potential of the vortex- free vector field j(x, t) Q. 

In linear approximation, when / — [Aijj^ + C(V'0)^]/2, we get from (4) and 

(5) 

Vj = T(Ai; - CA4,), (6) 

and j(k, w) — uj^{'k)S^(k,Lu)k/k'^ with ijj^(k) = — ir(yl + Cfc^). The condition 
j(k, a;)|]^^P = is guaranteed by Sip{'k,uj)\^^^ — , which follows from (1). We 
see that j consist of two transport modes: a dilatation one, which corresponds 
to an extension or a contraction of OP and a diffusion one, which corresponds 
to Ficken diffusion. If we suppose that the characteristic time of OP dilata- 
tion TA = 1/{TA) is much greater than the caracteristic time of OP diffusion 
Tc = l/{TCk^), i.e. fc2 >> 1/^2 (^ ^ ^y^jjj jg ^^^ correlation length of OP 
fluctuations), we get Pick's law in its classical form 

j = -VV^P (7) 



with the diffusion coefficient V = TC. 
To derive Pick's law in the form 



-AVm/, (8) 



where ^i — iJ.i{x, t) is a local chemical potential and A is a transport coefficient, 
it is necessary to introduce a local equilibrium assumption. Let us divide the 



system into cells of a small volume Vi — l"^ with I « 1/k. In each cell the 
diffusion mode of j(k, lj) dies away after the time r/ = P/{rC) much smaller 
than the relaxation time tc of the diffusion mode in the system as a whole. 
Then we can believe that all local thermodynamic variables and functions as 
OP -0, free energy Fi, entropy Si, pressure Pi, chemical potential /i; etc. are 
constants in the cell. It is the assumption of local equilibrium that makes 
possible to meaningfully define the local free energy Fi{x,t) that is the same 
function of the local thermodynamic variables as the equilibrium free energy is 
of the equilibrium thermodynamic parameters. That is why the fundamental 
differential relation dF; — —SidT — PidV + Vimdip is vahd locally. Replacing 
F{^p} = Jy^ <ix[AiP^ + C{V^)^]/2 by the local free energy F,('0) = ViAiP'^/2 and 
using formulas of traditional equilibrium thermodynamics we can calculate the 
local chemical potential iJ,i{x,t) = dFi{ip)/dtp = ViAilj{x,t) and find (8) with 
A = TC/A Q. If we are interested by self-organisation effects, we must consider 
the situation, where ^^fc^ « 1 and j is defined by (6), but not by (7) or by (8). 
We see that Pick's law in form (8) is itself a consequence of OP evolution 
equation, so it can not be utilised for its derivation. Pick's law in form (7) 
describes a special case of OP transport: diffusion without dilatation. Pick's law 
in form (8) is an approximative expression for diffusion flux, not only in the sense 
that it does not take into account the higher order gradient terms of a chemical 
potential /i; , but also in the sense that it is a law of the linear thermodynamics 
for which the local equilibrium assumption is compulsory. General expression 
for OP flux (5) has nothing to do with the Pick's law form 

j = -X\76F/S^P (9) 

criticized in [g and accepted in Q| as "motivated phenomcnologically" . The 
use of (8) with formally defined "local chemical potential" iJ,{x,t) — 5F/Sip has 
no thermodynamic foundation and is incorrect. This completely concerns the 
equation 

d^/dt = -XV^6F/di;. (10) 

An origin of the misunderstanding is, as has been note in 0], a false adop- 
tion of the Pick's law in form (9) as a fundamental general law, which can be 
used as a basis for derivation of the OP evolution equation of nonequilibrium 
thermodynamic systems with COP. 

Pick's law (8) fails completely in a non-linear dynamics. In a non-linear 
approximation, when /{V'} = [A'ijj^+BTp'^/2 + C(yip)^ + D{A'iP)^ + E{^\'^)^]/2 
equation (4) has the form 

Vj = T[Ai} + BV^ - CAtP + DAV - Ei^^AiP]. (11) 

We see that Vj includes now not only the dilatation mode and the diffusion 
mode, but also a cross term Eip'^Atp, being responsible for the dilatation - 
diffusion coupling. In this case we can not in principal separate dilatation and 
diffusion effects. 



Some words about the example given by Rutenberg. It does not show any 
inconsistence of my theory with the local consevation law (2), but brilliantly 
demonstrate incorrectness of (10). Let us, following W, take "a special initial 
condition" and "require that the dissipative dynamics be invariant under f/; —f 
~tjj and that F{tp} is minimized hy ip = ±1 everywhere except for a small 
sphere where ip = —1, the other of which has tjj = +1 and —1, respectively". 
For simplicity, let us omit the gradient term in Rutenberg's free energy funtional 
and take F{ip} = J dx{ip'^ — 1)^. Then the system initial state is evidently an 
equilibrium one for all three cases NCOP, GCOP or LOOP. Therefore we must 
not observe any changes of '0(x, t) in time. However, following ||l|, we get from 
(10) that "the spheres evolve". This is nonsense. This circumstance can be 
directly discovered from (10), if we note that by the condition V'^SF/S-ip = it 
defines a false equilibrium state, being different, in general case, from the real 
equilibrium one. The later must always minimize F{ip}, i.e. must be determined 
by the condition 5F/5ip — 0. 
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